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т.    я.    СБКЕРЖ-ЗБНЬКОВИЧ 

НЕКОТОРЫЕ  ЗАДАЧИ  ТЕОРИИ  РАСПРОСТРАНЕНИЯ 
ПРИЛИВНЫХ   ВОЛН  В  НЕОДНОРОДНОЙ  ЖИДКОСТИ 

1 

В  динамике  моря  мало  разработана  теория  распространения  прилив- 
ных волн  в  неоднородной  жидкости  при  наличии  вращения  бассейна. 

Простейший  учет  неоднородности  заключается  в  том,  что  жидкость 
представляют  состоящей  из  двух  или  нескольких  слоев  разной,  но  посто- 
янной в  каждом  слое  плотности.  Таким  образом,  задача  о  распростране- 
нии волн  в  неоднородной  жидкости  сводится  к  изучению  волн  на  свобод- 
ной поверхности  и  поверхности  раздела,  т.  е.  внутренних  волн.  Предла- 
гаемая работа  посвящена  некоторым  вопросам  теории  распространения 
длинных  волн  в  жидкости,  состоящей  из  двух   слоев  разной  плотности. 

I.  ВЫВОД  УРАВНЕНИЙ  ЗАДАЧИ  О  РАСПРОСТРАНЕНИИ 
ВОЛН  ВО  ВРАЩАЮЩЕМСЯ   БАССЕЙНЕ 

1.  Постановка  задачи  и  уравнения  движения 

Рассмотрим  бассейн,  наполненный  жидкостью,  состоящей  из  двух 
слоев  разной  плотности,  и  предположим,  что  бассейн  вращается  вокруг 
вертикальной  оси  с  постоянной  угловой  скоростью  со.  При  этом  будем 
считать,  что  бассейн  имеет  постоянную  глубину,  дно  бассейна  плоское, 
а  стенки  представляют  собой  цилиндрическую  поверхность  с  верти- 
кальными образующими.  Глубина  нижнего  слоя  жидкости,  наполняющей 
бассейн,  равна  /г^,  верхнего — /га",  плотность  нижнего  слоя  равна  р1,  верхне- 
го—рз- 

Предположим  теперь,  что  в  бассейне  распространяется  свободная 
волна.  Так  как  жидкость  двухслойная,  то  волны  образуются  как  на 
свободной  поверхности,  так  и  на  поверхности  раздела.  Исследуем  эти 
волны.  Для  этого  сначала  выведем  те  дифференциальные  уравнения, 
решая  которые  можно  будет  получить  вид  свободной  поверхности  и 
поверхности    раздела. 

Систему  координат  выберем  следующим  образом:  оси  Ох  и  Оу  распо- 
ложены в  плоскости,  параллельной  дну  бассейна,  и  участвуют  во  враще- 
нии бассейна;   ось  Ог  направлена  вертикально  вверх. 

При  выводе  уравнений  для  определения  возвышений  используются 
линейные  уравнения  Эйлера  в  подвижной  системе  координат.  Для  ниж- 
ней жидкости  они  имеют  вид: 


—  2сог;1  = 


^  +  2»«,= 


Р1      дх 

1    дрг 

Р1       дУ 


(1) 
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где  и^  и  1^1 —  компоненты  относительной  скорости  частиц  нижней 
жидкости,  направленные  по  осям  Ох  я  Оу;  рх  —  давление  в  нижнем 
слое.  Ускорение  силы  тяжести  входит  в  рх,  а  центробежное  ускорение  мы  в 
уравнения  не  ввели,  потому  что  в  дальнейшем  будем  отсчитывать  возвы- 
шение от  параболоида  враш;ения  как  формы  относительного  равновесия. 
В  третьем  уравнении  движения  считаем  вертикальное  ускорение 
малым  и  пренебрегаем  им,  т.  е.  рассматриваем  длинные  волны.  Из  этого 
уравнения,  как  известно,  получается  закон  гидростатики  для  опреде- 
ления давления: 

Р1  =  891  {К  +  Сх  —  21)  +  ^р2  (/«2  +  Са  —  Сх)  +  Ро,  (2) 

где  Сх  —  возвышение  поверхности  раздела  над  своим  средним  уровнем, 
который  является  параболоидом  вращения;  Сз  —  возвышение  свободной 
поверхности  над  своим  средним  уровнем,  также  являюш;имся  параболо- 
идом вращения,  2^  —  вертикальная  координата  рассматриваемой  точки, 
отсчитываемая  от  дна;  р^  —  атмосферное  давление;  §  —  ускорение  силы 
тяжести.  Глубины  /г^  и  /гд  мы  считаем    постоянными. 

Для  верхней  жидкости  уравнения  движения  будут  иметь  следующий 
вид: 


I 


-^  —  2сог;2  — ^ 

оЬ  ^  р2     ах 


(3) 


Здесь  проекции  скоростей  и  давление  обозначаются  теми  же  буквами, 
только  с  индексами  2  снизу.  Давление  также  определяется  по  закону 
гидростатики 

VI  =  8Н  (^2  +  ^2  —  22)  Н-  /?о, 

где  22 — вертикальная  координата  рассматриваемой  частицы,  отсчитываемая 
от  среднего  положения  поверхности  раздела. 

Подсчитывая  изменения  со  временем  количества  жидкости,  которое 
находится  в  элементарном  объеме,  получаем  уравнение  неразрывности 
для  нижней  жидкости: 

Ку     Ы  \  дх    '^    ду  )  ^  ^ 

и  для  верхней  жидкости: 

^^(С.-У  =  -(^  +  ^).  (5) 


ду 

Чтобы  получить  уравнения  нашей  задачи,  мы  поступим  следующим 
образом.  Рассмотрим  сначала  нижнюю  жидкость.  Из  уравнений  (1)  и  (4) 
мы  можем  получить  одно  уравнение  для  определения  возвышения  поверх- 
ности раздела  Сх-  Для  этого  продифференцируем  уравнение  (4)  по  времени: 


Лх    д1^    ~        [дх  \  д1  )  '^  ду  \  д1  }\ 


И  воспользуемся    при    преобразовании    правой    части  этого    соотношения 
уравнениями  (1).  Тогда  будем  иметь 


.    дК,  1    (д^Р1     ,    д^р,  \    ,    р^   /  аых  дь,  \  .^. 
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Но  ИЗ  (2)  получаем 

где  о  =  -^  ,  а  =  1  —  о,  А  —  оператор  Лапласа. 


Из  (1)  И  (4)  имеем: 

д  (  дщ 
д1  V  ду 

дх  ) 

2ы    3^1 

/11       д1 

Интегрируя  ПО  1,  находим: 

дщ 

9^1 

2а)  ^ 

ду 

дх     ~ 

'  Л1  ^^- 

(8) 

Добавочную  произвольную  функцию  интеграции,  зависящую  только  от 
X  и  у,  мы  отбрасываем,  потому  что,  как  можно  показать,  она  определяет 
установившееся   движение,    которое  в  данном    случае  нас  не  интересует. 

Подставляя  в  (6)  значения  7г(-^+^)и  (-^  —  -1^)^3(7)  и  (8), 
получаем 


дКх 


9,2  +4со2С1  =  С1^А(аС1  +  8С2),  (9) 

где  с1  =  §к^.  Уравнение  (9)  и  является  одним  из  основных  уравнений 
задачи.  Второе  основное  уравнение  получится  из  уравнений  (3)  и  (5), 
если  проделать  с  ними  аналогичные  преобразования.  Оно  имеет  вид: 

|^^  +  4соК2  =  с?Л(аС1  +  8С2)  +  с^АС2,  (Ю) 

где  Со  =  ё^2- 

Граничные  условия  для  интеграции  уравнений  (9)  и  (10)  получаются 
из  требования  равенства  нулю  скорости,  нормальной  к  стенкам  бассейна. 
Таким  образом,  имеем  для  нижней  жидкости: 

4^  +  2»§-  =  0.  (И) 

где  частные  производные  по  п  и  /  берутся  соответственно  по  направле- 
нию нормали  и  по  касательной  к  поверхности;  для  верхней  жидкости 
имеем: 

2.  Получение  окончательных  уравнений  задачи 

Таким  образом,  мы  теперь  имеем  для  Сх  и  Сг  систему  уравнений  (9) 
и  (10)  с  граничными  условиями  (И)  и  (12).  Для  решения  этой  системы 
удобнее  ее  преобразовать,  положив 

С1  +  5С2  =  2,  (13) 

где  5  есть  пока  неопределенный  постоянный  множитель.  Умножим  урав- 
нение (10)  на  5  и  сложим  с  (9).  В  результате  получим: 

^  +  40)22  =  ас?  (1  +8)М—  [ас?52  +  (ас?  -  Ьс\  -  с|)  5  -  Ьс1]  М^.       (14) 
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Множитель  5  определяем  так,  чтобы  удовлетворялось  соотношение 

асх^з  +  (ас1  —  к?  —  с1)8  —  Ьс1  =  0,  (15) 

т.  е.  чтобы  в  (14)  пропал  член  с  АСг-  Подробное  исследование  уравнения 
(15)  показывает,  что  его  корни  «1  и  ^2  действительны  и,  кроме  того, 
один  из    них,  например    корень    ^1,    всегда  положительный,  а  другой  «а, 

всегда  отрицательный,  но  |б'2|<1.  При  этом  предполагается,  что  о=  — 

Р1 

удовлетворяет    неравенству    О  <<  о  <^  1 ,  а  отношение  -^  изменяется    от    О 

"1 

до  оо.  Подставив  теперь  корень  5  —  ^1  в  уравнение  (14),  получим: 

-^  +  4со%  =  а?А^1;  (16) 


положив  в  том  же  уравнении  5  =  ^з,  находим: 

^^  +  46)2^2  =  «1А^2,  (17) 

где 

а^  =  ас?(1 +51),     а!  =  ас?  (1+52).  (18) 

Так  как  а1^0  и  а^'^О,  то  уравнения  (16)  и  (17)  являются  диффе- 
ренциальными уравнениями  в  частных  производных  гиперболического 
типа. 

Для  граничных  условий  получаем: 

И 

Таким  образом,  для  каждой  из  функций    Е^    и  ^2  имеем  отдельные  гра- 
ничные условия. 

Теперь  из  (13)  определим  1^  и  Сз-  Полагая  5  ==  ^1  и  5  =  ^з  и  вычитая 
одно  из  другого,  будем    иметь: 


^1 


^2^1 *1^2 

«2  ^1 


(19) 


^2  =  4^-  (20) 

51  ^2 

Следовательно,  задача  для  жидкости,  состоящей  из  двух  слоев  разной 
плотности,  сведена  к  двум  отдельным  задачам  для  однослойной  жидкости, 
так  как  каждая  из  функций  ^1  и  2>2  является  решением  отдельного  урав- 
нения распространения  волн  в  бассейне  постоянной  глубины.  По  найден- 
ным решениям  ^1  и  ^з  формулы  (19)  и  (20)  позволяют  найти  решение 
первоначальной  системы  уравнений,  определяющей  волны  во  вращаю- 
щемся бассейне  с  двуслойной  жидкостью.  Этот  метод  сведения  к  одно- 
родной жидкости  впервые  был  разработан  Л.  Н.  Сретенским  [И  для 
случая    двуслойной    жидкости,    покрывающей    сферу. 
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II.  НАХОЖДЕНИЕ  ЧАСТНЫХ  РЕШЕНИЙ  В  ЗАДАЧЕ  О  ВОЛНАХ, 
РАСПРОСТРАНЯЮЩИХСЯ  В  БЕСКОНЕЧНО  ДЛИННОМ 
ВРАЩАЮЩЕМСЯ  КАНАЛЕ  ПОСТОЯННОЙ  ГЛУБИНЫ 
С  ДВУМЯ  СЛОЯМИ  ЖИДКОСТИ 

Напомним  коротко  известные  результаты,  касающиеся  распростра- 
нения волн  во  вращающемся  канале  с  однородной  жидкостью. 

Пусть  ось  Ох  проведена  вдоль  одной  из  стенок  канала,  ось  Оу  идет 
перпендикулярно  оси  Ох  в  направлении  к  другой  стенке.  Рассмотрим 
частное  решение  уравнения  (16)  в  виде: 

^1  =  /1Ысо8(А;1а;— аО,  (21) 

дающее  при  фиксированном  у  прогрессивную  волну,  распространяющуюся 
вдоль  оси  Ох.  Здесь  к^ —  обратная  длина  волны,  а  —  частота.  Функция  Д  (г/) 
Определяется  из  дифференциального  уравнения 

/:(у)  =  Га?+-^^^^)Л(^).  (22) 


Обозначим 

^1  л 1 —  =  ^1  (23) 

«1 

и  будем  считать,  что  д'1^0.  Тогда,  интегрируя  уравнение  (22),  получим 

д  (г/)  =  Лгб^-У  +  В-^е-^^У.  (24) 

Произвольные  постоянные  интеграции  А^  и  ^1  определим  из  граничных 
условий,  которые  в  данном  случае  имеют  вид: 

а/;(0)  +  2б)Л1Л(0)  =  0  (25) 

и 

а/;(6)  +  2<оА:1А(6)  =  0.  (26) 

Подставив  Д  {у)    и    Д  {у)  из  (24)  в  граничные    условия    (25)    и    (26), 
получим 

(ад1  +  2(оА;1)  Л +  (20)^1 -0^0^1  =  0,  | 

е^.ь  (а^^  +  2о>А;1)  А-,  -\-  е-^'^  (2а>А;1  —  а^^)  Вг  =  0.]  \^  '> 

Чтобы  система  (27)  имела  ненулевое  решение,  определитель  ее  должен 
равняться  нулю.  Отсюда  находим 


(а^^  +  2б)А;1)  (2соА;1  —  ад^)  (е-^'^»  —  е^'Ь)  =  0. 
соотноше 
1)  а 


Это  соотношение  будет  выполняться  в  следующих  случаях: 


3)  е^^ь  =  е-9'Ь. 

Этому    последнему    случаю    отвечают    волны    Пуанкаре,    которые    мы 
в  данной  работе  изучать  не  будем. 


т.     я.     СЕКЕРЖ-ЗЕНЬКОВИЧ 


1 .  Рассмотрим    случай  а  = ^ . 

Тогда 

а2д1  =  46)4?.  (28) 

Подставляя  д1  из  (23)  в  (28),   имеем: 


(а2_4со2)А;1  = 
Отсюда,  если  о^фАа^,  получим 


А?  =  ^.  (29) 


Подставляя  Ах  из  (29)  в  (23),    имеем: 

^1  =  -1-  .  (30) 

Соотношения  (29)  и  (30)  выполняются  при  следующих  допущениях: 
А. 

7  а  2(0 

Тогда  из  (24),  положив  А^  —  С^,  8^  =  С^,  имеем 

Для    определения    констант  Сх    и   С^,  подставив   значения  А^,  В^^  к^ 
и  ^1  в  первое  уравнение  (27),   получим 

Следовательно,  Се  =  О,    С^  —  произвольно. 
Из  (21)  имеем  частное   решение: 

^1  =  С^е    «>  "^  С08  [^  —  оЛ  .  (31) 

Б. 

,  а  2со 

Положив  теперь  ^1  =  Сз,  Вх  —  С^,  получим,  так    же   как   и  в  преды- 
дущем случае,  частное  решение  вида 

2,[  :^  с;е  «'  "^  С08  ( -^  +  ^^)  .  •  (32) 

2.  Рассмотрим  случай    а  = ^ .  Таким  же   путем,  как  и  в  случае 

1,  мы  найдем  следующие  частные  решения: 

2сй  2ы 

21  =  С^е"  ^ ""  соз  [^  -  а  А  и    ^Г  =  С'ге^^\оз  [^  +  о{^  .  (33) 

Решения  уравнения  (17)  получаются  аналогично  и  имеют  вид: 

^2  =  С^е    ^ "  соз  [^  -  01^  ;      г^  =  С'^^^     соз  {—  +  а^^  .  (34) 
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Переходим  теперь  к  задаче  настоящего  раздела,  т.  е.  к  рассмотре- 
нию ВОЛН  во  вращающемся  канале,  содержащем  два  слоя  жидкости.  Пусть 
глубина  верхнего  слоя  жидкости  к^,  нижнего  кх  и  ширина  канала  Ь. 
Ось  Ох  проведем  вдоль  одной  из  стенок  канала,  ось  Оу  направим  пер- 
пендикулярно оси  Ох     в  направлении  к  другой  стенке. 

Выше  было  показано,  что  возвышения  поверхности  раздела  Сх  и 
свободной  поверхности  Са  определяются  формулами  (19)  и  (20)  через  ^^ 
и  2<2у  являющиеся  частными  решениями  уравнений  (16)  и  (17). 
Воспользовавшись  этими  формулами,  мы  получим  выражения  для  Сх  и  Сг- 
Мы  будем  брать  только  суммы  двух  волн  типа  Кельвина  с  одинаковыми 

27Г 

периодами  —  . 

Рассмотрим  следующие  случаи: 

1.  Возьмем  частное  решение  уравнения  (16)  и  частное  решение  урав- 
нения (17),  соответствующие    волнам    типа    Кельвина    противоположного 

направления,  с  разными  амплитудами  и  длинами  (/^1= и  Ьг  =  "V"; 

Подставив  их  в  формулы  (19)  и    (20),  получим 

2м  ,  .  2<о 

у.Схе     °»    С08  (  -^ <^И  —  51^2  е  "•     соз  ( -^  -Ь  о  Л 

—  ■—V  [ах  \  ,   —У  I  ^^     ,       Л 

^  ^1^     °'     соз(— -с.^]-С,в°»      сов(— +  а^) 


'2 


«2 


2.  Взяв  два  частных  решения  уравнения  (16)  для  случаев  суммы  двух 
встречных  волн  типа  Кельвина,  имеющих  разные  амплитуды,  но  одина- 
ковые длины,   получим: 


2со  ,  ^  2ы 


^^    ^  .,  [С,е     "'  '  соз(^  -  <^^)  +  С[еа.  "  .р,  (^  +  ы)] 
—  —  V  /ах  \  I   7ГУ        I  ^^    ,      Л 


*1  *'2 

3.  Берем  два    частных    решения    уравнения    (17)    опять    для  того  же 
случая,  что  и  в  п.  2.  Подставив  их  в  формулы  (19)  и   (20),  найдем: 

2й>  ,  ,  2ш 


^    _   -1  [Сге     °'  '  соз  (-^  -  ог)  4-  С;е  "■  '^  соз  (^  +  а/)] 

52  — «1 

2со 

оЛ  л-  С'^е'^г  ^  соз  [^  +  оП 


52  — «1 
2со  ,  ,  2со 

V         [ах 

С^е     °»     соз  ( — 

г      \  ^2 

^•2  


*1  *2 


4.  Взяв  частное  решение  уравнения  (16)  и  частное  решение  уравне- 
ния (17),  соответствующие  волнам  типа  Кельвина  одинакового  направле- 
ния, но  с  разными  амплитудами  и  длинами,  и  подставив  их  в  формулы 
(19)  и  (20),  получим: 
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2и  ,  .  2<о 


~7ГУ         /стж            \          /-.— г-У         А  ста:  \ 

52^16     «'       соз(-^^ стП  — «Ае     «2      сов! стП 


«2  51 

2(0  2ы 


С^в     «.      соз^— -ст^)-С2е     «^     соз(— -ст^) 


С2  = 

Отметим,  что  распространение  парных  комбинаций  волн  Кельвина 
и  Пуанкаре  на  поверхности  однородной  жидкости  во  вращающемся  ка- 
нале   изучал    Фишо    [2]. 

III.  ИССЛЕДОВАНИЕ  ПОЛУЧЕННЫХ  ЧАСТНЫХ  РЕШЕНИЙ. 

ИЗУЧЕНИЕ  ПОЛОЖЕНИЯ  АМФИДРОМИЧЕСКИХ  ТОЧЕК 

И  ПОСТРОЕНИЕ  КОТИДАЛЬНЫХ  КАРТ 

Рассмотрим  подробно  только  первый  из  указанных  в  конце  раздела 
II  случаев,  т.  е.  волны,  распространяющиеся  в  противоположных  направ- 
лениях и  имеющие  разные  амплитуды  и  длины.  Для  остальных  случаев 
ограничимся  приведением  результатов. 

1.  Волны  на  поверхности  раздела 

Исследуем  свободную  поверхность  и  поверхность  раздела  в  окрест- 
ности некоторых  особых  точек  —  около  амфидромических  точек  и  около 
точек   с   горизонтальными   касательными   плоскостями. 

А.  Исследование  области  около  амфидромических  точек 
Возвышение  в  этом  случае  дается  формулой: 

2со               ,  ^  2(л  ,  '         . 

8^3уе     "'      соз стП  —  зуС^  е"»      соз    -—  +  ^Ч 

г      ^  \^^1 /  ^ \    «2  /  /35) 

"1  5,-81 

Для  определения  положения  амфидромических  точек  преобразуем 
выражение  (35)  следующим  образом: 

2й>  2ы  . 

/  чг  (     /-1     ~1ГУ  ож  уо'~^  сгжЛ  ,, 

(^2  —  ^1)^1  =  (^2^16    "'     С08 ^^Сгб"»    С08  — 1соза?  + 

,  2со  2ы^  . 

+  иС^е    ^  %1п  -^  -Ь  1С;е  «»  ^  зхп  ^)  8Ш  а^  (36) 

В  амфидромической  точке  Сх  должно  быть  тождественно  относительно 
I  равно  нулю,  т.  е. 

С1  =  0.  (37) 

Из  (36)  видно,  что  (37)  выполняется,  если 

2со                                              2со 
^    —  -—  У          ста;              п'  'ТГ'^          °^         п 
5оС,е      «'     СОЗ 51^26"'      СОЗ =  О, 

2м  2ы  ох  {  \Р^) 

^ У    •      стж     ,        /->'  ~^  У    •     л 

52^16      «'      81П \-  8уСф°-^      ЗШ  аз    =  <^- 
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Так  как  С^  и  С^  отличны  от  нуля,  то  определитель    системы  должен 
быть  равен  нулю,  т.    е. 


2о  2<о 

—  —-  У            ох           —-У  ах 

е     **'      С08  —  е  "•  СОЙ  — 

2й>  2<о 

—  Г—  У    .       ох           7Г~У  •      ох 
е      "'       81П е "»  8111  — 

«1  «2 


=  0. 


Отсюда  получаем: 


X  = 


жМ      а1а2 


о        «1  +  «2    ' 


(39) 


где  Л^  —  любое  целое  число.  Подставляя  в  первое    уравнение  (38)  значе- 
ние X  из  (39),  имеем: 


2.,  («_1±Л1)  у 


^     С08 
*2  С\  ^\  4"  Яг 

2  С08 


(40) 


Но 


а\  +  ^2 


С08 

■кNа_ 

а^Лг^г 

соя 

■кНа^ 

(-1Л 


^1  "Г  Я^ 

поэтому  г/  определяется  из  (40)  в  виде: 


2/=  Т-Г^^-т1п[(-1)^^'^ 


(41) 


Из  (41)  видно,  что  если  С^  и  С^  одного  знака,  то  Л^  надо  брать  не- 
четным, если  же  они  разных  знаков,  то  N  должно  быть  четное,  так  как 
^1  и  ^2  разных  знаков.  Таким  образом,  мы  видим,  что  на  поверхности 
раздела  имеется  бесконечное  множество  амфидромических  точек. 

Определим  теперь  вид  котидальных  линий  около  амфидромических 
точек.  Обозначим 


2<о  2со 

—У           ОХ  '    ТГ  ^  ^^ 

^2^16      "»      С08  — ЗхС^е^*       С08  — 


«1 


«2 


=  л, 


2<й 


2йз 


*2  —  *1 


=  в,. 


Тогда  из  (35)  имеем  С1  =  А^  С08  о1  -|-  ^1 81п  о1  или 

г^  =  асо8(а^  —  а). 


если  положить 


А^  =  а  С08  а,     3^  =  а  81п  а,     а  =  агс  1§ 


^^1 


(42) 


(43) 


(44) 


Котидальными  линиями  называются  линии  одинаковой  фазы  прилива. 
Следовательно,  их  уравнение,  в  силу  (43),  имеет  вид 


о1  —  а  (ж,  у)  =  С. 


(45) 
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Фиксировав  I  и  меняя  константу  С,  получим  семейство  котидальных 
линий  для  данного  момента  времени,  т.  е.  котидальную  карту  этого- 
момента  времени. 

п 

На  основании  (45)  и  введя  обозначение  1^  =  1 ,  имеем 


X  {х,  у)  =  а1  —  С  =  а({ 1  =  а^^. 


Отсюда  в  силу  (44)  и  отбрасывая  индекс  у  ^1,  получаем  уравнение 
котидальных  линий  в  виде 

|1  =  Ь§а^  (46) 

Подставив  сюда  значения  А^  и  В^  из  (42),  находим: 

5„С,  (?    °''    зт ^8,  С„е"2    51П  — 

^ ^ =  ^ё  ^«.  (^7). 

—  —  У        ста;  ^'     Г~  У        ох 

5„С,  е    °^^     соз •  8,  С,  е  "«     соз  — 

2    1  й!  •^    -=  Со 

После  разложения  в  степенные  ряды  функций,  стоящих  в  числителе 
и  знаменателе,  в  окрестности  амфидромической  точки  с  координатами 
(^0)  ^о))  ^  также  на  основании  (38)  уравнение  (46)  можно  записать  в  виде- 


где 


Ь{х-хо)+с{у-~у,)  =  0,  (48) 

О  =  а  I  — —  е  со8  — ^^  +  — —  е    '       соз  — !^  )  + 

\   «1  а-,  йп  Со    / 


«1  «1  «2  «2 

^             2ш  /  2ы 

,          .           ,|^2^"1  ~~^^''  •      ОХо            ^1^^2  "^  ^''     .       ОЖо 

^               2<о  /  2С0 


V         «1  «1      '      йа  «2  у 

+  2со1да^и-1е    о.       со8-^+-^е^У»со8-^   . 

\    "1  "1  "2  ^2  / 


Уравнение  (48)  показывает,  что  около  амфидромических  точек  коти- 
дальные  линии  представляют  собой  прямые,  проходящие  через  амфидро- 
мическую  точку.  Определим  теперь  направление  вращения  котидальных 
линий.  Из  (38)  имеем: 

^     ~  ^  ^         ста:  ^'    ~^  У  ах 

З^Сув  соз —  =51626  С08 , 

,  Ях  Яг 

2ы  2(0 

^    ~  ^'аГ  ^    •     ста;  ^'    ~^  ^    .     ох 

з^Слв  зт  —  =  —  5]Сое  зш  —  . 


Отсюда 


2й> 

у 


СОЗ =  7*1^2<?  "'     »  (49) 


2со 
V 


СОЗ =  7^2^!^        "'        »  (50) 
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8111 ==  —  ЬзХов 

а,  ^    2 


Оз 


Отсюда 


Возведя  в  квадрат  (49)  и  (51)  и  сложив,  получаем: 

1 


7^  +  0^  = 


4Сй 


Следовательно , 


Т  = 


1                  ,      .                1 
— С08  ф,     О  = 

2<о  Т'  5 


'А' "' 


^^;«"^ 


81Пф. 


Из  (50)  и  (52)  аналогично  находим: 


(51) 
(52) 


(53) 


а, 


С08ф,      8  = 


8Ш  ф. 


З1С1  ^1^1 

Подставив  (53)  и  (54)  в  (49),  (50),  (51)  и  (52),  получим 


СОЗ =С03ф:      81П 

«1  '  «1 


зшф;  СОЗ 


созф;  з1п  —  =  зшф. 

'  «2 


(54) 


(55) 


В  силу  (55)  для  Ъ  и  с  будем  иметь  выражения: 

г  Л2С1 ^^''      ,     ^1<^2      ~^^'\ 

Ь  =  а[  -^-^  е      '       +  -^^-^  е    '         соз  Ф  — 

2со  „'        2м 


\     01 


V     <^1 


_  о  ,  2со 


Разделив  первое  из  этих  равенств  на  второе,  находим 

Ь       с(1— 1д(тМдф) 


с  2(0  (1дф  4-  1^0^)  * 

Вычисляем  производную  по  и 

(Ц  ^ ^  ^8^  Ф  +  1 

\С  /1  2С|>  СО! 


<о. 


005%^      {Ь§'\>  +  1^о1)^ 

Следовательно,  —  убывает  с  ростом  I.  Из  (48)  имеем: 

Ь 


У  — У  о 


(х  —  Хо). 


(56) 
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Вращение  котидальных  линий  происходит  против  часовой  стрелки, 
так  как  —  6/с  —  тангенс  угла  наклона  котидальной  линии  уменьшается 
по  абсолютной  величине. 

Найдем  линии  равной  амплитуды  около  амфидромических  точек.  Эти 
линии  определяются  уравнением 


Ч 


4  +  5?  -  с\ 

Подставив  в  (57)  значения  А^  и  В^,  находим: 


(57) 


[     /-,           а,                ста;             /-,'      а,    ^           ож  \2    , 
5,616  С08 «Хоб      '       С08 + 

2со  2со 


(58) 


После  разложения  в  степенные    ряды   выражений,  стоящих  в  скобках 
и  используя  (55)  и  (38),  запишем  (58)  в  виде: 


2"  „/        2С0 

Уо  5,  С„      -I-  Уо 


2"  /        2Ь) 


2«о  _,/        2со 

-:7-  Уо  5,  С. 


\  «1 


+  '-11?.^")в1п.>.(,-,,)|  =  С^ 


Положив  ^е      "'  ''°  +  — е  ''^  ''°  =  Р,   из  (59)  имеем: 

а,  '      ао 


(59) 


Фиг.  1 


а2  (:г  -  :Го)2  +  4«2  •  (г/  -  г/о)^  =  С  * ,  (60) 

где     С*  —  новая     константа,     равная 

-— .     Уравнение   (60)    представляет 

собой  уравнение  семейства  подобных 
эллипсов.  Следовательно,  около 
амфидромических  точек  семейство 
линий  равных  амплитуд  есть  семей- 
ство подобных  и  подобно  располо- 
женных эллипсов. 

На  фиг.  1      изображена  амфидро-    I 
мическая     точка     Л.    От    нее  расхо-    ' 


дятся  котидальные  линии  —  сплош- 
ные. Линии  равных  амплитуд  изображены  пунктиром.  Стрелки  ука- 
зывают направление  вращения  котидальных  линий  около  амфидромиче- 
ской  точки. 
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В.  Исследование  области  около  точек  с  горизонтальными 
касательными    плоскостями 

Найдем  на    поверхности    раздела  точки  с    горизонтальными  касатель- 


д-С.^ 


д^г 


ными  плоскостями;  для  этих  точек  -^—  =  У  и  -^ 

ох  ду 

Из  (36)  ВЫЧИСЛИМ   -^  : 


0. 


2со 


+(^.Ав 


2со 
1—  У 


С08 8лЬф  С05 —    8т  а^. 

«1     '     Ог     ^  аг 


Приравняем    -^  нулю.  Это  равенство    будет    выполняться  при  любом 


I,  если 


201 

■::-  V 


з^Слб  8т ^Хоб    "     8т  —  =  о, 


о         ^ 


а,                ох            о         .1'      а,    "           ох 
С08 ^Хоб  С08 

0102  «2 


0. 


(61) 


Система  (61)  с  неизвестными  Сх  и  Сг  имеет  не  нулевые  решения,  если 
определитель  ее  равен  нулю,  т.  е. 


I 


2со 
Г-  У 


8,>е 


«1 


Зов 


81П 


С08 


8ле 


8,е 


8т 


С08- 


-0. 


Отсюда  получаем: 


X  = 


ТГ^У         «102 


_^      ,  (62) 

<  где  уУ  —  любое  целое  число.  Значение    у   определяем  из  первого  уравне- 
ния (61). 

!У-  .    П     ,1пГ(-1)^+^^^41.  (63) 

*         Отсюда  видно,  что    если  С^  и  6*2    одного  знака,  то  N  надо  брать  чет- 
1    ным,  если  же  они  разных  знаков,   то  N  должно  быть  нечетное. 

^  Вычислив  теперь    производную  -^    и  приравняв  ее  нулю,  мы  получим 

для  X  и  г/  те  же  значения  (62)  и  (63).  Это  значит,  что  существуют  такие 

точки,  в  которых  обе  производные  -^  и  -^  равны  нулю,    т.    е.    точки  с 

горизонтальными  касательными  плоскостями.  Котидальные  линии  около 
этих  точек  представляют  собой  семейство  гипербол  с  асимптотами,  па- 
раллельными осям  координат.  Их  уравнение  имеет  вид: 


-  2соа  [-1-  С1е      "'   '°  +  8,8,С,С',  [\-^\)^ 
+  -^С,е^'^■      \ху=^Е. 


(^-ё)у* 


+ 
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Линии  равных  амплитуд  около  точек  с  горизонтальными  касательными 
плоскостями  являются  семейством  гипербол  с  асимптотами,  наклоненными 
к  осям  координат.  Они  определяются   уравнением: 


где 


—  аЦх  —  Жо)2  +  4со2  {у  —  г/о)2  =  С*, 


з^С^е 


у« 


+  чС'ге  "' 


"2 


Уо 


И  С  —  произвольная  постоянная. 

На  фиг.  2  изображены  котидальные  ли- 
нии (сплошные)  и  линии  равных  амплитуд 
(пунктирные)  около  точки  с  горизонтальными 
касательными  плоскостями  (точка  Б). 

В  результате  проведенного    нами    иссле- 
дования можно  построить  для  данного    слу- 
чая   волн    на    поверхности  раздела,    например,   следующую  котидальную 
карту,  которая  представлена  на  фиг.  3  б.  На  ней  изображены  котидаль- 
ные линии  (сплошные)  и  линии    равных    амплитуд    (пунктирные). 


Фиг.  2 


2.   Волны  на  свободной  новерхности 

В  данном  случае  подробного  исследования  приводить  не  будем,  а 
дадим  только  результаты,  так  как  необходимый  разбор  делается  анало- 
гично только  что  проведенному.  В  рассматриваемом  случае  возвышение 
дается  формулой 


^2  = 


2<о  , 

^    — ;;"  ^       /  «уа; 

Схв       "1       соз ст^ 

\а1 


«1  «2 


(64) 


В  результате  исследования  мы  получили   следуюш;ие   координаты  ам- 
фидромических  точек: 

У  = 


X  = 


а^а^, 


а     01  +  «2  ' 


„    ,      ,      ,  1п— ;5- ; 

2о>(а1-Ьа2)         С 


N  —  четное,  если  С-^  и  С^  одного  знака. 

Котидальные  линии  около  амфидромических  точек  являются  прямыми, 
враш;ающимися  около    амфидромических    точек    против  часовой  стрелки. 

Линии  равной  амплитуды  представляют  собой  семейство  подобных 
и  подобно  расположенных  эллипсов. 

Координаты  точек  с  горизонтальными  касательными  плоскостями  на 
свободной  поверхности: 

^  ^  теЛ^     а^а^  ^  а-^а^  ^^  <^1    Дд 

ст    «1  +  аг  '  ^         2(0  (Ох  +  «г)         С'^    «1 


НЕКОТОРЫЕ    ЗАДАЧИ    ТЕОРИИ    РАСПРОСТРАНЕНИЯ    ПРИЛИВНЫХ    ВОЛН        17 

Котидальные  линии  около  этих  точек  являются  семейством  гипербол 
с  асимптотами,  совпадающими  с  осями  координат.  Линии  равной  ампли- 
туды представляют  собой    семейство    гипербол   с    асимптотами,  поверну- 


Фиг.  3 

а  —  котидальная  карта  для  свободной  поверхности;  б  —  то  же  для  поверхности 

раздела 

тыми  относительно  координатных  осей.  Котидальная  карта  для  рассматри- 
ваемого случая  изображена  на  фиг.  3,  а. 

3.  Сравнение  возвышений  свободной  новерхности 
и  поверхности  раздела 

Представим  возвышения  поверхности  раздела  и  свободной  поверхности 
в  виде 


:.^  =  Уа1-\-в1соз{о1-{-о^1), 

Сз  =  Уа1  +  В1  С08  {аЬ  +  аз). 


Пусть  отношение  амплитуды  возвышения  свободной  поверхности  к 
амплитуде  возвышения  поверхности  раздела  в  одной  и  той  же  точке 
{х,  у)  равно  некоторому  числу  к,  т.  е.        ' 


А1  +  В1 


(65) 


Подставив  в  (65)  значения  А^  и  ^1  из  (42)  и  аналогичные  выражения 
для  ^2  и  ^2.  которые  можно  получить  из  (64)  и  проделав  простые  пре- 
образования, находим 


_2<о(— 4--^)у  „„        2»  ( \-—)у 

=  2С,С',{к\8, - 1)  соз  <:х (^  +  ^)  . 


(66) 


Для  исследования    этого    уравнения    введем   новые   переменные  с  и  тг]  по 
формулам: 


«=«(^+^> '"=24^+4)^- 


2     Труды  МГИ,  т.  VIII 
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Тогда  (66)  примет  вид 

С1  (1  _  1с^81)  е-''  +  С'1  (1  —  к^81)  е"  =  1С^С'^  {к\8^  -  1)  со8  е.       (67) 

Решаем  это    уравнение    относительно    т];    здесь    возможны  следующие 
случаи: 

1.  к  =  0.  Из  (65)  видно,  что  тогда  А2  =  0  и  ^2  =  0.  Эти  соотношения 
определяют  амфидромические  точки  на  свободной  поверхности. 

2.  А;  =  оо.  Из  (65)  следует,  что  ^^  =  О   и   ^х  =  О,  т.  е.  получаем  ам- 
фидромические точки  на  поверхности  раздела. 

3.  А;<' — ,  т.  е.  к  очень  мало,  так  как  ^^  велико. 
Решаем  уравнение  (67)  в  этом  случае.    Обозначим 

{1-к'81)С1^^Яе-\ 


{1  —  к^81)С'1  =  ^Яе\ 


где 


7?  =  2]/  (1  -  к^81)  (1-^81)  С  А,  а  =  у1п- 


п>2, 
^2 


причем  а>0,  если  — ^>-1.  Тогда  уравнение    (67)   представится    в  виде 


—  Л  сЬ  (т^  -\-  (X.)  =  т^  С08  8, 


(68) 


где 


т'^ 


2{1-к\8,)СА. 


Уравнение  (68)  решаем  графически.  Для  этого  строим  кривые  —  /?  сЬ  (тг)  +  а) 
и  т^соз|,  откладывая  аргументы  ^  и  т]  на  оси  абсцисс,  а  значения  функ- 
ций—  на     оси  ординат,    как 
показано  на  фиг.  4. 

Берем  точку  А'  кривой 
т^со5^1  с  абсциссой  ^1,  рав- 
ной ОА,  и  из  этой  точки  А' 
проводим  прямую,  параллель- 
ную оси  абсцисс  до  пересе- 
■"  чения  с  кривой  —  Л  сЬ  (т;  -}-а) . 
Получим  две  точки  пере- 
сечения В  ж  С.  Затем 
от  точки  А  на  прямой,  па- 
раллельной оси  ординат, 
откладываем  отрезки,  равные 
О'В  и  О'С.  Получим  две 
точки  К  ж  Ь,  изображаю- 
ш;ие  корни  уравнения  (68) 
для  ^  =  ^1.  Затем,  взяв  точки  кривой  т^соз^  с  абсциссами  ^2>  ^з>  •  •  •  >  мы 
проводим  те  же  построения  и  в  результате  получаем  замкнутую 
кривую  в  виде  овала.  На  этой  кривой  отношение  амплитуд  возвышений 
равно  к.  Внутри  овала  Д<<А;Л1.  Вне  овала  А2'>кА1,  т.  е.  внутри 
овала  амплитуда  поверхности  раздела  будет  в  1/к  раз  больше  ампли- 
туды свободной  поверхности,  вне  —  в  1/к  раз  меньше  (здесь  к  очень 
мало). 


^ 

к 

^  г 

^ 

\/ 

А/ 

/' 

;\ 

\ 

г 

— 

^\ 

л 

А 

0> 

с^^ 

л} 

к 

1 

Фиг.  4 
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4.  А;^  >  -5- .  т.  е.  к  велико.  В  этом  случае  получаем  овалы  около  ам- 


фидромических    точек    на    поверхности    раздела.    Внутри    овалов    имеем 
Л^'^  кА^.  Вне  овалов  А2<^кА1. 


5.  к^  =  ^. 


6.  А^  =  -^ 


Для  случаев   5    и  6  проводим  построения,  аналогичные    предыдущим 
(фиг.  5).  Получаем  11-образные  кривые. 


Фиг.  5 


7.  Общий    случай,    т.    е.    имеем    1 — к^81^0,  1 — к^81<С0  и 
—  1  +  к^8182  <С  О'  Обозначив 

С1  (1  -  к' 81)  =  \  Яе-\  С'^  (1  -  А^^?)  =  -  \Ке\ 

где  /?  ^  О  и  а  >  О,    имеем 


где 


8Ь  (т)  +  а)  =  т^  С08  ^, 
т2  =  2  (1  —  кН^8^)СгС'.,. 


Решая,  как  и  выше,  это  уравнение  графически,  получим  волнистую 
кривую,  изображающую  г^  как  функцию  \. 

Эта  кривая  делит  плоскость  на  две  части.  В  той  части,  которая  рас- 
положена выше  этой  кривой,  амплитуда  свободной  поверхности  больше 
амплитуды  поверхности  раздела  в  к  раз,  в  другой  части  амплитуда 
свободной  поверхности  в  к  раз  меньше  амплитуды  поверхности  райдела. 

Давая  к  различные  значения,  будем  иметь  семейство  таких  кривых. 
При  /с  =  1  получим  кривую,  выше  которой  амплитуда  поверхности 
раздела  больше  амплитуды  свободной  поверхности,  а  ниже  кривой  — 
наоборот.    Вдоль    этой  кривой  амплитуды  будут  равны. 

Окончательные  результаты  этих  исследований  изображены  на  фиг.  6. 

Следует  отметить,  что  так  как  ширина  канала  может  быть  произволь- 
ной и  стенки  его  могут  быть  расположены  произвольно  (но  параллельно 
прямой   ММ)   по   отношению    к    рисунку,    изображенному  на  фиг.  6,  то 
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ВОЗМОЖНЫ  следующие   случаи   в  зависимости   от    расположения   берегов 
канала: 

1.  В  канале  нет  областей,  где  к'^1.  Это  будет  в  том  случае,  если  оба 
берега  расположены  выше  максимумов  кривой,  отвечающей  А;  =  1. 


Фиг.  6 

А,    Ао    Аз,  А,  —  амфидромнческие    точки  на  свободной  поверхности;  В,,  'Вз,  Вз  —  то  же  на 
оов'ерхности  раздела.    Заштрихована  область,    где  амплитуда   возвышения   поверхности  раз- 
дела больше  амплитуды  возвышения  свободной  поверхности 


2.  Оба   берега   расположены    ниже   минимумов   кривой,    отвечающей 
к  ^?  1,т.  е.  в  областях,  где  А;>1.  В  этом  случае  в  канале  не  будет  областей, 


Фиг.  7 

а  —  котидальная  карта  для  свободной  поверхности;  б  —  то  же  для   поверхности  раздела 

Где  амплитуда  возвышения  поверхности  раздела  будет  больше  амплитуды 
возвышения  свободной  поверхности. 

3.  Один  из  берегов  пересекает  кривую,  отвечающую  А;  =  1.  В  канале 
будут  области,  в  которых  амплитуда  поверхности  раздела  больше  ам- 
плитуды   свободной   поверхности,    а   также   будут   области,    в    которых 


I 
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амплитуда  поверхности  раздела   меньше   амплитуды   свободной   поверх- 
ности. 

В   заключение   настоящего  раздела  приводим  котидальные  карты  для 
двух  других  случаев.  Карты  на   фиг.  7  относятся   к  случаю,  когда  рас- 


Фиг.  8 

а  —  котидальная  карта  для  свободной  поверхности;  б  —  то   же  для  поверхности  раздела 

сматривается  сумма  двух  встречных  волн,  имеющих  разные  амплитуды, 
но  одинаковые  длины.  Карты  на  фиг.  8  относятся  к  случаю,  когда  рассмат- 
ривается сумма  двух  волн  одинакового  направления,  имеющих  разные 
амплитуды   и  длины. 

Из  рассмотрения  карт  видно,  что  карты  для  поверхности  раздела  и 
свободной  поверхности  во  всех  случаях  оказались  качественно  подобными. 


11У.  ЗАДАЧА  ОБ  ОТРАЖЕНИИ  И  ПРОХОЖДЕНИИ  ВОЛН, 
РАСПРОСТРАНЯЮЩИХСЯ  В  КАНАЛЕ  ПЕРЕМЕННОЙ  ГЛУБИНЫ, 
К010РЫЙ  СОДЕРЖИТ  ДВА  СЛОЯ  ЖИДКОСТИ  РАЗНОЙ  ПЛОТНОСТИ 

1.  Распространение  волны 
из  глубокой  части  канала  в  мелкую 

Рассмотрим  бесконечно  длинный  канал,  состоящий  из  двух  частей 
разной,  но  постоянной  глубины.  Жидкость  в  его  глубокой  части  состоит 
из  двух  слоев  разной  плотности;  при  этом  нижний  слой  доходит  до  уровня 
дна  второй  части  канала,  а  верхний  слой  переходит  в  мелкую  часть.  Плот- 
ность верхнего  слоя  рг,  нижнего  р1.  Глубина  верхнего  слоя  Лг,  нижнего  кх. 

Считаем,  что  канал  не  вращается  и  решаем  одноразмерную  задачу. 
Положение  оси  Ох  указано  на  фиг.  9.  Ось  Оу  направлена  перпендикулярно 
стенкам  канала. 

Предположим,  что  из  отрицательной  бесконечности  по  поверхности 
канала  распространяется  свободная  волна  в  положительном  направ- 
лении оси  Ох,  как  указано  на  фиг.  9. 

Дойдя  до  плоскости  раздела  х  =  0,  часть  волны  отразится,  а  часть 
пройдет  далее,  в  мелкую  половину.  Так  как  в  канале  имеется  два  слоя 
жидкости,  то  на  поверхности  раздела  также  образуются  волны,  которые, 
дойдя  до  плоскости  раздела,  отразятся  от  нее.  Задача  состоит  в  том,  чтобы 
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исследовать  волны  как  на  свободной  поверхности,  так  и  на  поверхности 
раздела. 


Фиг.  9 


Для   определения   возвышений  используем   выведенные   в   разделе-  1'. 
уравнения  (9)  и  (10);  однако  в  них  теперь  пропадут  члены  с  производными 
по  у,  так  как  решается  одноразмерная  задача  и,  кроме  того,  о)=  О,    по- 
тому что  канал  не  вращается. 

Таким  образом,  уравнения  для  определения  возвышений  имеют  вид: 


д1^ 
д1^ 


дх^ 


.2^ 

^1  дх^ 


2  дх^ 


О, 


(69) 


где  Сх  —  возвышение  поверхности  раздела;  Сг  —  возвышение  свободной 
поверхности  в  глубокой  части  канаЛа;  Сз  —  возвьппение  свободной  поверх- 
ности в  мелкой  части  канала;  8  =  -^,  с^  =  {1—Ь)§к^^,    с^  =  дк^. 

Третье  уравнение  системы  (69)  —  одноразмерное  волновое  уравнение. 
Решаем  систему  (69).  Так  как  в  первые  два  уравнения  входят  только 
Сх  и  Сз,  а  в  третье  только  Сз,  то  решаем  отдельно  систему  первьк  двух 
уравнений. 

Она  имеет  частное  решение,  которое  дает  прогрессивные  волны  вида: 


С.  =  ^^(*-Т).  ^1=В^('-Т> 


(70) 


где /^и ]  —  данная  произвольная  функция,  отвечаюш;ая  волне,  рас- 
пространяющейся в  положительном  направлении  оси  Ох.  Коэффициенты 
А  ж  В  связаны  так,  чтобы  удовлетворялись  первые  два  уравнения  сис- 
темы (69).  Подставив  Сх  и  Сг  из  (70)  в  эти  уравнения  и  приведя  к  общему 
знаменателю,  получим: 


§\ЬА  —  (с2  —  с2)  5  =  О, 

(с2  —  с2  —  дк^Ь)  А  —  с\В  =  0. 


(71) 


Система  (71)  имеет  отличное  от  нуля  решение,  если  определитель 
системы  равен  нулю.  Раскрывая  этот  определитель  и  делая  приведение 
подобных  членов,  получим  биквадратное   уравнение    для   определения  с: 


с' -8  {кг  +  к,)  с2  -\- §Ц1  -  Ь)кгк,  =  0. 


(72) 
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Уравнение  имеет  четыре  корня  {±^с;  гЬ^'),  причем  два  корня  попарно 
равны  по  абсолютной  величине,  но  противоположны  по  знаку;  далее 
считаем,  что  |  с  |  >  |  с'  |.  Следовательно,  система  (69)  имеет  общее  решение, 
состоящее  из  суммы  частных  решений,  которые  представляют  прогрес- 
сивные волны: 

С1  =  В/'((-^)  +  В'/-((+^)  +  Я/'(*+^),    [  (73) 


'  -КРИ-^);  ) 


-3 


здесь  Сз  —  частное  решение  третьего   уравнения. 

Считаем,  что  из  бесконечности  распространяется  данная  волна  с  амп- 
литудой А,  а  остальные  коэффициенты  неизвестны.  При  этом  коэффи- 
циенты имеют  следующий  смысл:  А'  ж  В'  —  амплитуды  отраженных  волн, 
идущих  с  большей  скоростью  с  по  свободной  поверхности  и  поверхности 
раздела;  V  и  Е  —  амплитуды  отраженных  волн,  идущих  с  меньшей  ско- 
ростью с'  по  свободной  поверхности  и  по  поверхности  раздела;  К  —  амп- 
литуда проходящей  волны,  идущей  со  скоростью  Сд. 

Рассмотрим  теперь  граничные  условия  задачи  и  уравнения  для  оп- 
ределения неизвестных  коэффициентов.  При  х  =  О  должны  выполняться 
следующие  граничные  условия: 

1.  Условие  непрерывности  возвышения,  т.  е. 

С2  =  Сз,  (74) 

или,  подставляя  в  (74)  значения  Сд  и  Сз  из  (73),  будем  иметь: 

А  +  А'-]-В  =  К.  (75) 

2.  Условие  непрерывности  горизонтальных  скоростей,  т.  е. 

«2  =  Щ  (76) 

для  любого  момента  времени. 

Из  уравнений  горизонтального  движения  жидкости  скорости  были 
выражены  через  возвышения: 

ди2  _  о^^       ^  ^  ^^ 

д1   ~        ^  дх   '     д1  ^  дх  '      . 

Но  так  как  скорости  равны  в  любой  момент  времени,  то 


ИЛИ 

дх  дх 


(77) 


Подставляя  в  (77)  значения  Сг  и  Сз  из  (73)  и  приводя  к  общему  зна- 
менателю, будем  иметь: 


ее' 


с'  (Л'  -А)  +  сП  =  --^К.  (78) 
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3.  Горизонтальная  скорость  течения  нижней  жидкости  равна  нулю 
в  любой  момент  времени,  т.  е. 

щ  =  0; 
следовательно, 

^  =  0.  (79) 

Выразив  щ  через  Сх  и  Сг,  представим  (79)  в  виде 

(1-«)§  +  2§  =  0  (80) 

Подставляя  значения  Сх  и  Сг  из  (73)  в  (80),  получим 

(1  —  8)  [с'  {В'  —  В)'^сЕ]+Ь  [с'  {А'  —  А)-\-  сВ]  =  0.  (81) 

Кроме  уравнений  (75),  (78)  и  (81),  имеем  еш;е  три  уравнения  для 
определения  неизвестных  коэффициентов,  а  именно:  первое  уравнение 
системы  (71)  и  аналогично  еще  два  уравнения: 

^^18^'  —  (с2  —  с1)  В'  =  О,  (82) 

§к^ЬВ  —  {с^  —  с1)Е  =  0.  (83) 

Таким  образом,  получим  шесть  уравнений  с  шестью  неизвестными, 
а  именно,  первое  уравнение  (71)  и  уравнения  (75),  (78),  (81),  (82),  (83). 

Для  решения  системы  уравнений  (75),  (78),  (81),  (82),  (83)  и  первого 
уравнения  (71)  прежде  всего  выразим  В,  В'  и  ^  через  А,  А'  та.  В  из  урав- 
нений (82),  (83)  и  первого  уравнения  (71): 

В  =  а^А,   В'  =  а^А',1Е  =  а^В,  •           (84) 

где 

8             1                           8             1  /огч 

А         1  —  8  _21_1                       1  —  8     с^  .                                   ^      ' 

Подставив  теперь  В,  В'  ж  Е  ий  (84)  в  уравнения  (75),  (78)  и  (81), 
получим: 

А'  ^В  —  К  =  —А, 

с'А'  +  сВ-^  —  К  =  с'А,  (86) 

с'Ьт_А'  +  сЬ^П  =  с'Ь^_А, 
где 

^1  =  «1  —  Ьа-^  -{-  Ь,  Ь^  =  (12  —  8^2  4"  8.  (87) 

Решая  систему  (86),  будем  иметь 

Л'  =  "^  "^  [{с  +  сз)  Ъ,  +  {с-  с^)  Ь^] ,  (88) 

где 

А  =  '-^  [(С2  +  с')  Ь,-{с  +  с^)  Ь,] ;  (89) 

0^^^2^Ь,,  (90) 

К=^-^сс'{Ь,-Ь,).  (91) 
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Подставляя  значение  А  из  (89)  в  (88),  (90),  (91)  и  деля  числитель    и 
знаменатель  на  Ъ^с^,  находим  из  (88),  (90),  (91)  и  (84): 


_ц \с, ^д        д^ 


2 

с 

С2 

с' 

+  1- 

( 

С 

Ч 

-) 

/ 

2     1 

\ 

— 

Ь2 

-) 

Со  V    С2  У    01  / 


(92) 


Нам  надо  вычислить  амплитуды  Л\  В,  В,  В',  Е  и  К  для  разных  от- 
ношений глубин,  которые  являются  переменными  величинами.  Для  этого 
можно  было  бы  выразить  корни  биквадратного  уравнения  (72)  через 
отношение  Л^/Ла,  подставить  их  в  выражения  для  амплитуд  (92)  и  затем 
вести  расчет,  но  это  очень  сложно.  Лучше  ввести  безразмерный  вспомо- 
гательный параметр  6,  через  который  просто  выражаются  корни  уравне- 
ния (72)  и  коэффициенты  (92).  Затем  значения  параметра  Ь  надо  выра- 
зить через  отношение  кх/Ь,^-  Тогда  окончательно  получим  требуемую 
зависимость.  Итак,  поступаем  следующим  образом.  Из  биквадратного 
уравнения  (72)  имеем: 

сМс^-С2)  +  с?(с^-ТС^)  =  0,  (93) 

гдет  =  т=У 
Отсюда 


^1 

-^= 1^-  (94) 


^2 


Теперь  положим 
Тогда  из  (94) 


4=^-  (9^> 


'^^        '-^  (96) 


г-5 


Воспользовавшись  формулой  для  суммы   корней   квадратного   уравнения, 
получим  из  (93),  используя  (96): 

^-1  =  5.  (97) 


Найдем  также 

с'2         1  —  5 


(98) 


Теперь  из  формул  (85)  и  (87)  можно  выразить  а^,  «а,  Ъ^  и  бд  через  %. 
Получим  следуюпще  значения: 

л       ь  Т  — 1        и         Г  — ^  Л 

и    .._.  С(т-1)         Ь,  5(у-1)  ( 
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Подставив  найденные  значения  в  (92),  имеем: 
У     у  — Е  "^  /I/ 


(Г^  +  1)(у-С) 


/^ 


^    ,    ^  П(т-1) 

5  (КТ-1)(т-?) 


-4, 


^)-  — 2 


У 


1-^ 

г-5 


1^5  (у -1) 


^, 


(^^5-1)(у-5) 


Б  =  (1— е)Л,  Л'  =  (1-^)Л',  Е 


5(т-1) 


а:=2 


(5-1)(Т-5) 


/I 


5—1       ' 
_  Л. 


+  1 


1^ит-1) 


у-5    '    ^        (Г5-1)(у-5) 
По  этим  формулам  будем  определять  коэффициенты    как   функции    %. 

1-? 


Из  формулы 


т-5 


можно  установить,  что   надо    брать   следую- 


щие значения  Ъ 

Кроме  того,  между  Н-^/  к^  и  ^  существует  соотношение: 

Л1  _     1  —  5        1 

/га    "■'^  Т-^    *     ^    • 
0<5<1 


Т  аб лица  1 


с| 

/11 

В 

А' 

В' 

в 

Е 

к 

^        с, 

Л. 

А 

А 

А 

А 

А 

А 

0 

оо 

1 

1 

1 

—0,27 

10,65 

1,73 

0,2 

4,02 

0,8 

0,43 

0,35 

—0,15 

7,19 

1,29 

0,4 

1,52 

0,6 

0,27 

0,16 

—0,09 

5,77 

1,22 

0,6 

0,68 

0,4 

0,15 

0,06 

—0,04 

4.7 

1,15 

0,8 

0,26 

0,2 

0,06 

0,01 

—0,01 

2,93 

1,05 

1 

0 

0 

0 

0 

0 

0 

1 

у<5<оо 


40 

42 

44 

46 

48 

50 

100 

150 

200 

250 

оо 


—39 

0,73 

19,52 

41 

0,78 

9,77 

43 

0,81 

6,52 

—45 

0,33 

4,89 

47 

0,84 

3,92 

—49 

0,85 

0,66 

—99 

0,96 

0,36 

—149 

0,98 

0,25 

—199 

0,99 

0,19 

—249 

0,99 

0 

— оо 

1 

—28,32 

—32,16 

—34,76 

—37,24 

—39,59 

—41,89 

—94,67 

—145,54 

—196,10 

—245,44 

—  оо 


О 

-0,32 
-0,48 

0,42 
-0,53 
-0,57 
-0,85 
-0,91 
-0,94 
-0,95 
-1 


О 
—0,30 
—0,38 
—0,42 
—0,44 
—0,45 
—0,34 
—0,24 
—0,18 
—0,15 

О 


1,73 
1.47 
1,39 
1,38 
1,39 
1,28 
1,10 
1,07 
1,05 
1,04 
1 
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Если  построить  график  этой  функции,  то  можно  видеть,  что  двум  раз- 
личным значениям  $  соответствует  одно  и  то  же  отношение  к^!!!^,  при- 
чем большему  корню  с  отвечают  значения  ^,  изменяющиеся  в  интервале 
0^^^-^!,  а  меньшему  корню  с'  соответствуют  значения  ^,  изменяюга;иеся 

В  интервале  '^<^^<^оо.  Это    видно  из    формулы  -|- =  |,  Поэтому   нужно 

различать  два  случая.  В  первом  случае  данная  волна  идет  со  скоростью, 
отвечающей  большему  корню    биквадратного  уравнения,  во    втором   слу- 


^ — I — 1—] I I \ I I I I I 1__1 I \ I I I    I  I 1_^ 1-4, >_-2 


Фиг.  10 


чае  данная  волна  идет  со  скоростью,  отвечающей  меньшему  корню  этого 
уравнения. 

Были  проведены  вычисления  для  случая  -у  =  40,  что  соответствует 
отношению  плотностей  8  =  39/40.  Результаты  вычислений  приведены 
в  табл.  1.  Эта  таблица  состоит  из  двух  частей.  Первая  часть  составлена 
для  значений  ^  в  интервале  от  О  до  1 ;  вторая  часть  —  для  значений  ^ 
в  интервале  от  -^  до  оо.  Эти  значения  8  и  соответствующих  Лх/Т^а  приве- 
дены в  первых  столбцах;  в  остальных  столбцах  даны  значения  отношений 
амплитуд  волн  к  амплитуде  А  заданной  волны. 

На  основании  табл.  1  построены  графики  (фиг.  10  и  И).  При  пост- 
роении графика  фиг.  10  использована  табл.  1    для  ^  от  О  до  1.    На   оси 
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X  откладывается  отношение  глубин  /гх/Лг-  По  оси  ординат  откладываются 
величины  амплитуд  в  долях  А.  График  фиг.  10  построен  для  случая, 
когда  заданная  волна  идет  с  большей  скоростью  с.  Когда  глубина  ниж- 
него слоя  мала  по  сравнению  с  глубиной  верхнего,  т.  е.  отношение  /гх/Л^ 
близко  к  нулю,  амплитуды  всех  волн  очень  малы,  меньше  амплитуды 
данной  волны  А,  за  исключением  К  амплитуды  волны,  проходящей  на 
мелкую  часть  канала.  Величина  К  будет  немного  больше  А. 


/<-^1.73Д 


д^^-39Д 


Фиг.  И 

При  увеличении  к-^  /  к^  очень  сильно  растет  Е  —  амплитуда  отражен- 
ной волны  на  поверхности  раздела,  идущей  с  меньшей  скоростью  с'. 
Когда  /г^/Лд-^оо,  т.  е.  глубина  нижнего  слоя  значительно  больше  глу- 
бины верхнего,  имеем  ^-^10,6^.  Амплитуда  /С  при  увеличении  к^/Н^ 
растет  и  в  пределе  стремится  к  1,73  А.  Амплитуды  остальных  волн  при 
этом  растут  медленно,  все  время  оставаясь  меньше  А. 

График,  изображенный  на  фиг.  11,  относится  к  случаю,  когда  данная 
волна  идет  со  скоростью  с',  равной  меньшему  корню  уравнения  (72)  и  фор- 
мулы (73)  принимают  вид: 

1,  =  АР{1—^)-^А'р{1  +  ^)  +  Ор[1  +  ^-), 
1^^  ВР{1  -^)  +  В'Р[1  •^■^)  +  ЕР  {г  -^■^)  , 

(.,  =  Ке[1-- 
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Для  построения  графика  взята  вторая  часть  таблицы  для  значений 
$  от  40  до  оо.  На  ОСИ  абсцисс  откладываем  отношение  глубин  к^/к^. 
На  оси  ординат  откладываем  амплитуды  в  долях  А. 

Когда  глубина  нижнего  слоя  мала  по  сравнению  с  глубиной  верх- 
него слоя,  т.  е.  отношение  кх/ к^  близко  к  нулю,  то  очень  большие 
значения  будут  иметь  абсолютные  величины  \В\  и  \В'  \  амплитуд  пря- 
мых и  отраженных  волн  на  поверхности  раздела,  идущих  с  меньшей 
скоростью  с'. 

При  увеличении  кх/к^  эти  амплитуды  быстро  уменьшаются  и  при 
Л^/Лд— ^о©  будут  иметь  следующие  значения: 

В-^  —  2>^А,  Б'->  — 28,32Л. 

Остальные  волны  при  всех  отношениях  к^/к^  будут  иметь  незначи- 
тельные амплитуды  по  сравнению  о,  В  ж  В'  (амплитудами  волн  на  по- 
верхности раздела,  идущими  с  меньшей  скоростью). 


2.  Распространение  волны 
из  мелкой  части  канала  в  глубокую 

Условия  задачи  в  этом  случае  точно  такие  же,  как  в  п.  1.  Отличие 
заключается  в  том,  что  данная  волна  идет  с  мелкой  части  канала  на 
глубокую. 

Для  определения  возвышений  мы  используем  уравнения  (69). 

В  данном  случае  эта  система  уравнений  имеет  общее  решение,  кото- 
рое состоит  из  суммы  частных  решений,  представляющих  прогрессивные 
волны: 


Считаем,  что  данная  волна  распространяется  с  амплитудой  А,  осталь- 
ные коэффициенты  неизвестны. 

Граничные  условия  точно  такие  же,  как  в  п.  1.  Таким  же  способом, 
как  в  предыдущем  случае,  мы  получаем  уравнения  для  определения 
неизвестных  коэффициентов  и  так  же  их    решаем.    Затем    введя,    как    и 

в    первом    случае,    вспомогательный     безразмерный    параметр     5  =  -|- , 
получим  следующие  выражения  для  коэффициентов: 


П(т-1) 


^  =    Г-. — ;г- чгтт:. — 7Т ^' 


/^ 


5   ^1  П(г-1) 


(П  +  1)(Т-?) 
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В  = 


КТ(Т-1) 


(Г-?)(?-1)[|/1- 

Г      у  — ^ 


~1  ^  ^      у;^(у-1) 

5  (П-1)(у__^)] 


^, 


^1  =  2 


1/1-^   ^1 ^^(Т-1) 


^, 


т-^ 


(П-1)(у-5) 


^)  =  (1-1)5,  ^),  =  x^^5,. 


По  этим  формулам  проводились  вычисления,  результаты  которых 
приведены  в  табл.  2  для  значений  ^  в  интервале  0-<^<]1. 

При  значениях  ^  в  интервале  от  '^(^  до  оо  величины  амплитуд  будут  та- 
кими же,  как  и  для  промежутка  от  О  до  1,  так  как  в  данном  случае  перво- 
начальная волна  всегда  идет  с  одной  постоянной  скоростью  сг.  Так  же 
как  и  ранее,  считаем  ^  =  40. 

В  остальных  столбцах  табл.  2  даны  значения  отношений  амплитуд  волн 
к  амплитуде  А   заданной  волны. 

Таблица   2 


с1 

Л, 

А, 

В 

в, 

в 

Ох 

^-    с. 

/12 

А 

А 

А 

А 

А 

0 

схэ 

—0,73 

0 

0,27 

0 

-10,65 

0,2 

4,02 

—0,29 

0,57 

0,15 

0,45 

—  7,16 

0,4 

1,52 

—0,18 

0,76 

0,09 

0,46 

—  5,88 

0,6 

0,68 

—0,10 

0,84 

0,04 

0,34 

—  4,30 

0,8 

0,26 

—0,05 

0,94 

0,016 

0,19 

—  3,12 

1 

0 

0 

1 

0 

0 

0 

На  основании  табл.  2  построен  график,  изображенный  на  фиг.  12. 
В  этом  случае  данная  волна  амплитуды  А  идет  со  скоростью  с^-  На  оси 

абсцисс  откладываем  отношение  глубин  -^.  На   оси  ординат  откладываем 

амплитуды  в  долях  А. 

Когда  отношение  /г^  /  к^^  близко  к  нулю,  все  волны  имеют  малые  амп- 
литуды, меньшие,  чем  амплитуды  данной  волны  А. 

При  увеличении  к^  /  к^  очень  быстро  растет  \01\  —  амплитуда  волны 
на  поверхности  раздела,  идущей  с  меньшей  скоростью.  Когда  к^/ к^-^оо, 
1)1-»  — 10,65^;  ^5  — амплитуда  волны,  проходящей  с  большей  ско- 
ростью с  на  свободной  поверхности,  при  увеличении  к^/к^  уменьшается, 
стремясь  к  нулю. 

Амплитуда  В  волны  на  поверхности  раздела,  идущей  с  большей  ско- 
ростью, при  увеличении  Лх/Лд  сначала  растет,  а  затем  убывает  и  в  пре- 
деле стремится  к  нулю. 

Амплитуды  остальных  волн  при  росте  ^1/^2  увеличиваются  не  на- 
много, оставаясь  все  время  меньше  А. 
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В  заключение  рассмотрим  случай,  когда  жидкость  нижнего  слоя  глу- 
бокой части  не  доходит  до  дна  мелкой  части  канала  (фиг.  13).  Сохраняя 
обозначения  предыдущих  разделов,  будем  считать,  что  нижняя  жидкость 
с  плотностью  рх  не  доходит  до  уровня  дна  мелкой  части  канала.  Тогда, 
очевидно,  задача  сведется  к  решению  системы  уравнений  (69)  со  следую- 
щими граничными  условиями  при  гс  =  0: 


Л,^^-т.В5А 


С2    = 


Фиг.  12 

1.  Условие  непрерывности  возвышения,  т.   е. 

2.  Условие   равенства   расходов:    кгич.  =  /гзмз. 

3.  Горизонтальная  скорость 
течения  нижней  жидкости  равна 
нулю,    т.  е.  Й1  =  0. 

Решение  системы  (69)  при 
этих  граничных  условиях  может 
быть  проведено  тем  же  способом, 
что  и  в  предыдущих  разделах. 
Здесь  также    можно    рассмотреть 

два    случая:   во-первых,  заданная  ^^^-  ^^ 

волна  идет  с  глубокой    части  ка- 
нала на  мелкую  и,  во-вторых, —  с  мелкой  части  на    глубокую. 


Ш///////////////}. 


32  Т.  Я,  СЕКЕРЖ-ЗЕНЬКОВИЧ 

В  заключение  отметим,  что  результаты,  полученные  в  этом  разделе, 
можно  применить,  например,  к  исследованию  приливов  в  устьях  фиордов 
и  рек,  впадающих  в  моря  или  океаны. 
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